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ABSTRACT 

According to recent results of iHo & Heink3 ( |2009|) and iHeinke & Hoi ( |2010|) , the Cas- 
siopeia A supernova remnant contains a young 330 yr old) neutron star (NS) which has 
carbon atmosphere and shows noticeable decline of the effective surface temperature. We re- 
port a new (November 2010) Chandra observation which confirms the previously reported 
decline rate. The decline is naturally explained if neutrons have recently become superfluid 
(in triplet-state) in the NS core, producing a splash of neutrino emission due to Cooper pair 
formation (CPF) process that currently accelerates the cooling. This scenario puts stringent 
constraints on poorly known properties of NS cores: on density dependence of the temperature 
Tcnip) for the onset of neutron superfluidity [T^nip) should have a wide peak with maximum 
(7 - 9) X 10^ K], on the reduction factor q of CPF process by collective effects in super- 
fluid matter {q > 0.4), and on the intensity of neutrino emission before the onset of neutron 
superfluidity (30-100 times weaker than the standard modified Urea process). This is serious 
evidence for nucleon superfluidity in NS cores that comes from observations of cooling NSs. 
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1 INTRODUCTION 



NS cores contain superdense 
are still uncertain (see, e.i 



It is well known that 

matter whose properties are s t ill uncertain (see , e.g .. 
iHaensel. Potekhin & YakovlevI |2007| : iLattimer & Prakasd 120071) 
One can explore thes e properties by studying the cooling of 
isolated NSs (see, e.g.. |Pethicklll992l: lYakovlev & Pethick|[2004l ; 



IPage. Geppert & WebeJl2006l : |Page et al.ll2009l , for review) 

We analyse observations of the NS in the supernova rem- 
nant Cassi opeia A (Cas A). The distance to the remnant is d = 
3.4+55-5 kpc teeedetalj|l995h . The Cas A age is reliably estimated 
as t X 330 ± 20 yr f rom observations of the remnant expan- 
sion ( lFesenetal.ll2006l) . The compact central source was identi- 
fied in first-light Chandra X- r ay obs er vations jTananbaum [T999I) 



and studied by Pavlov et al. I 1 I2OOOI) : IChakrabartv et al.l feOOlh : 
IPavlov & Lun3 l l2009l) but its nature has been uncertain. The fits of 
the observed X-ray spectrum with magnetized or non-magnetized 
hydrogen atmosphere models or with black-body spectrum re- 
vealed too small size of the emission region (could be hot spots 
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on NS surf ace although no pulsations have been observed, e.g., 
IPavlov & Luna..2009). 

Recenfly IHo & Heink3 l l2009h have shown that the observed 
spectrum is successfully fitted taking a carbon atmosphere model 
with a low magnetic field (B < 10" G). The gravitational mass of 
the object, as inferred from the fits, is M a 1.3 - 2Mq, circum- 
ferential radius R x 8-15 km, an d the non-redshifted e ffective 
surface temperature ~ 2x 10* K jYakovlev et al.ll2oTll) . These 
parameters indicate that the compact source is an NS with the car- 
bon atmosphere. It emits thermal radiation from the entire surface 
and has the surface temperature typical for an isolated NS. It is the 
you ngest in the fami l y of o bserved cooling NSs. 

lYakovlev et al.l ( 1201 ll) compared these observations with the 
NS cooling theory. The authors concluded that the Cas A NS has 
already reached the stage of internal thermal relaxation. It cools via 
neutrino emission from the stellar core; its neutrino luminosity is 
not very diff'ere nt from that provided b y the modified Urea p rocess. 

Following IHo & Heinkj ( |2009|) . iHeinke & Hoi ( 1201 Ol) anal- 
ysed Chandra observations of the Cas A NS during 10 years and 
found a steady decline of by about 4%. They interpret it as di- 
rect observation of Cas A NS cooling, the phenomenon which has 
never been observed before for any isolated NS. These results are 



© 2010 RAS 



2 P. S. Shternin et al. 



confirmed by new observations we report below. We interpret them 
as a manifestation of neutron superfluidity in the Cas A NS. 

When t his paper was nea rly completed we became aware of 
the paper by IPage et aL who proposed similar explanation 

of the Cas A NS observations. However the two papers are differ- 
ent in details, and can be regarded as complementary. In particular, 
we discuss the dependence of cooling curves on the poorly known 
efficiency of neutrino emission due to CPF process and the possi- 
bility to interpret observations of all cooling stars by one model of 
superdense matter. It is important that we report the new observa- 
tion. 



2 NEW CHANDRA OBSERVATIONS 

W e use the Chandra da ta on the Cas A NS, discussed and fitted 
bv lHeinke&Hol l llOlOl) . and add one new data point. Briefly, the 
analysed data include the Chandra ACIS-S observations of Cas A 
longer than 5 ks. In order to ensure that all considered data are 
directly comparable, we take only the 'directly comparable data' of 
iHeinke & Hoi (|2010|). W e exclude the suban'ay obser vation of Cas 
A ^Pavlov & Lunall2009l) . since the pileup properties ( lDavisll200lh 
of this spectrum differ from the others, and those observations in 
which the Cas A NS dithered over bad pixels (most of the 2004 
observations). 

The new data point is pr oduced from two ACIS-S observa- 
tions of Cas A ( |Patnaudeetal]|2010 : ObsIDs 10936, 13177) taken 
on 2010 October 31 and 2010 November 2 for 33 and 17 ks (re- 
spectively) and telemetered in GRADED mode (as done for pre- 
vious Cas A observations). We used CIAO 4.2 and CALDB 4.2.1 to 
reprocess the data and produce response functions, correcting for 
the time-dependent ACIS quantum efficiency degradation and gain 
changes, but not for the charge-transfer inefficiency (as this cannot 
be modeled with GRADED data). We extracted source spectra with 
a 4 pixel (2.37") radius reg ion, and background spectra from an an- 
nulus of 5 to 8 pixels, as in lHeinke&Hoi l l20I0l) . We combined the 
two new observations into one set of spectra and responses (with 
an effective date of 2010 November 1), and grouped the spectrum 
by 200 counts/bin. 

We fitted these spectra simultaneously, forcing the NS mass 
and radius, along with the distance an d A^h. to be the same as in 
the best fit of [lleinke & Ho ( 20i3) and lYakovlev et al.l ( l20Ilh . and 
finding the Icr errors on the surface temperature at each epoch. This 
fitting is designed to cleanly define the relative variation in the tem- 
perature, separating this question from the absolute uncertai nty in 
the temperature (described in detail in I Yakovlev et al. I l201lh . The 
fitted values of the non-redshifted effective surface t emperature Ts 
(Table 1) are slightly different from those reported bv lHeinke & Hoi 
( I2OIOI) for the 2004-2009 observations, but within the Icr errors. 
The key result is that the new observation confirms and extends the 
cooling trend seen in lHeinke&Hol ( l2010l) . 



3 SECOND TEMPERATURE DROP 

The observed surface temperature decline is too steep and can- 
not be described by the theory discusse d bvlYakovlev et al. Ilmil) 
(who did not analyse the decline itself). iHeinke & Hoi ( |20IO|) sug- 
gested that Cas A NS undergoes the last years of the internal 
crust-core relaxation accompanied by a pronounce d surface tem- 
perature drop. However, the theory predicts (e.g.. [ Lattimer et al.l 
1 1994 iGnedin. Yakovlev & Potekhi3l200lh I Yakovlev et alj|20 III) a 



Table 1. Carbon atmosphere s pectr al fits, using the best spectral fit (M, 
R, Nh) of iHeinke & Hoi j2O10l) and lYakovlev et al.l j201lh . with the addi- 
tion of 2010 data. Epoch dates are for the midpoints of the observations, or 
weighted midpoints of merged datasets. Temperature errors are Icr confi- 
dence for a single parameter 



Epoch 


Exposure 


logr. 


ObsID(s) 


(Year) 


ks 


K 




2000.08 


50.56 


f. a7co+0.0019 


114 


2002.10 


50.3 


6 3237+™'^ 


1952 


2004.11 


50.16 


6 3156+™'' 


5196 


2007.93 


50.35 


6.3108S|^ 


9117, 9773 


2009.84 


46.26 


6.3087;0:««|« 


10935, 12020 


2010.83 


49.49 


6.3O60S15 


10936, 13177 



shorter relaxation, lasting typically < 100 yr. We propose another 
interpretation based on the effects of supe rfluidity in NS cores, in 
line with the work o f iGusakov et"ai] ( l2004 . 

Neutrons, protons (and other baryons if present) in NS interi- 
ors can be superfluid (due to Cooper pairing). Free neutrons in the 
inner crust and protons in the core undergo Cooper pairing in spin- 
singlet state, while neutrons in the core can pair in spin-triplet state. 
Critical temperatures of superfl uidity onset 7'c(p ) are very model 
dependent (as reviewed, e.g., bv lPage et al.l2004l) . 

Our idea is that the initial crust-core relaxation in Cas A NS is 
over, leaving the star sufficiently warm. We assume further that not 
too strong triplet-state neutron superfluidity appeared in the core 
some time ago. It initiates a splash of neutrino emission due to CPF 
process producing the second drop that mimics the second (de- 
layed) thermal relaxation. The second drop is mainly regulated by 
(i) the density dependence = T^^ip) of the critical temperature 
for the neutron pairing in the NS core, (ii) the reduction factor q 
of CPF process by collective effects in superfluid matter and (iii) 
the neutrino luminosity prior to the onset of neutron superfluidity. 
These properties are very model dependent and poorly known but 
can be constrained from the Cas A NS observations. 

The CPF process was first predicted by 

iFlowers. Ruderman & Sutherland ( 1 1 9761) for singlet-state 
pairing neglecting collective superfluid effects. Collec- 
tive effects can suppress t he pr ocess (q < I) as was 
first noticed by [Leinson (2001) and calculated la ter b y 
Leinson & Perej 1I2OO6); S e drakian. Muther & Schuckl l l2007h 
Kolomeitsev & Voskresenskvn2008l) : Steiner & Reddy ^ 
Leinsod ( I2OI0I) : some discussion is also given by IPage 



2009); 



et 



1 200^ : the results are controversial. The main attention has 
been paid to CPF process due to singlet-state pairing of neutrons 
which has been found to be strongly suppressed (q <k 1). We are 
interested in triplet-state pairing in which case the suppression 
is thought to be less dramatic. If the latter CPF process were not 
affected by collective effects, then 24% of neutrino emissivity 
would go through the vector weak interaction channel, while the 
rest (76%) would go through the axial-vector channel. Collective 
effects suppress the vector channel almost completely but the 
axial-vector channel survives. Exact value of q is debatable (the 
lowest estimate q w 0.19 is given bv iLein son 2010). Instead of 
relying on any specific model, we consider ^ as a free parameter. 

We illustrate this idea by cooling simulations of NSs 
with nucleon cores. We take the Akmal-Pandharipande- 
Ravenhall (APR) equation o f state (EOS) in the core 
l lAkmal. Pandharipande & Ravenhal l, 1998j. By APR we mean the 
param etrization of APR results by iHeiselberg & Hiorth-Jenseiil 
i ll999l) - the version APR I proposed bv lCusakov et al.l ( l2005l) . The 
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p[io"gcm'] t[yr] 

Figure 1. (Color on line) Left: Models (a), (b) and (c) for critical temperature of triplet-state neutron pairing in NS core. Vertical dotted lines show central 
densities of NSs with M = 1 and 1 .65 M© (for APR EOS). Thin sohd line is the temperature profile in the 1 .65 Mq star of age 330 yr with neutron superfluidity 
(a). Right: Lines (a), (b) or (c) show cooling of the 1.65 Mq star with strong proton superfluidity and moderate neutron superfluidity (a), (b) or (c) in the core 
(q = 0.76) compared with observations of Cas A NS cooling. The dotted line is the same for the 1.75 Mq star and neutron superfluidity (a). The line (a) for 
M = 1.65 Mq is also shown in the inset over longer time span. Lines pSF and N in the inset are calculated for M = 1.65 Mq neglecting, respectively, neutron 
superfluidity and entire nucleon superfluidity in the core. 



maximum mass of a stable NS for this EOS is Mn,;,x = 1.929Mq; 
the powerful direct Urea process of neutrino emission is open in 
stars with M > M^u = L829Mo. Let the direct Urea process 
and even less efficient modified Urea process be either not 
ailowed or strongly suppressed in the Cas A NS. Otherwise the 
star would be too cold after the initial crust-core relaxation; we 
wouid be unable to significantfy speed up its coofing by the CPF 
process. To suppress Urea processes we assume the presence of 
strong proton superfluidity in the core, with criticaf temperature 
Tcpip) <; (2 - 3) X 10' K; exact vafues of T^p are unimportant. It 
occurs within a few days after the star formation. The proton CPF 
neutrino emission does not influence the cooling e ven if this emis- 
sion were not affected by collective effects (e.g., lYakovlev et al.l 
l200lh . As long as neutrons are non-superfluid, the neutrino 
emission is mainly generated in neutron-neutron bremsstrahlung 
process. It is wealc and leaves the Cas A NS sufficiently warm be- 
fore the second temperature drop. Unless the contrary is indicated, 
we consider NSs with ordinary ( non-accretted, non-magnetized) 
heat blanketing envelopes (e.g., IPotekhin. Chabrier & YakovlevI 
Il997l : IPotekhin et alJboOSi) and neglect superfluid efl'ects in the 
stellar crust (w hich weakly affect the cooling after the crust-core 
relaxation; e.g.. lYakovlev et al]|200lh . Note, that a splash of CPF 
neutrinos makes the star slightly non-isothermal, with the cooling 
slightly dependent on the thermal conductivity in the core, but the 
isothermal state is soon restored. 

Illustrative results are shown in Figs. [T] |2] and [3] The left 
panel of Fig. [T] gives three models of T^^ip). We do not rely on 
any specific theoretical model but consider three phenomenolog- 
ical curves (a), (b) and (c) (locating neutron superfluidity at pro- 
gressively lower densities). In our case the maximum of TcnCp) is 
strictly constrained to Ten max « (7 - 9) x 10** K. Higher or lower 
Tcnip) peaks would start the second temperature drop in the Cas A 
NS earlier or later than required by the observations (or even com- 
pletely wash out this drop in a very young or old star). The rcn(p) 
profiles should not be too narrow and the reduction factor q should 
not be too small (otherwise the second temperature drop would be 
weak). Other plots give theoretical cooling curves T^(t) (i.e., the 
redshifted surface temperature versus t). 

The right panel of Fig. [T] gives three cooling curves for the 



1.65 Mq star of nearly Cas A NS age with T^^ip) models (a)-(c) 
from the left panel and with q = 0.76 (axial vector channel is not 
suppressed). The curves are compared with the Cas A NS data. 
All three superfluidity models agree with these data. Taking su- 
perfluidity (a) and increasing the stellar mass to 1.75 Mq makes 
the CPF process more important and starts the second tempera- 
ture drop earlier (the dotted curve), in disagreement with the data. 
However, we could easily readjust (slightly decrease) T^^ip) and 
explain the Cas A NS data with the 1.75 Mq model. Therefore, we 
can successfully fit the data for a range of NS masses and take the 
1.65 Mq star as an example. The inset shows the cooling curve (a) 
(M = 1.65 Mo) over longer time scale. The second temperature 
drop at / ~ 250 yr is clearly pronounced here. Also, we show the 
cooling of non-superfluid star (curve N) and the star with proton su- 
perfluidity alone (curve pSF). Without the CPF process, the slopes 
of both curves are much smaller than required. 

In Fig. |2] we test two neutron superfluidity models, (a) and 
(c), against observations of other isol ated NSs. Observational data 
are taken from refere nces cited in lYakovlev e t al. (200^ and 
iKaminker et"aLl ( |2009|) with the exception of PSR J0007 -1 -7303 . 
The data on the latte r source are taken from lCaraveo et aL 
and | Lin et al. Note that in similar fig. 5 of lYakovlev et al.l 

l l201lh source labels 12, 13 and 14 should refer to RX J1856-3754, 
Geminga and PSR B 1055-52, respectively. 

We make a reasonable assumption that all NSs have different 
masses but the same physics of matter in their cores. We present 
cooling curves (solid lines) for a number of masses (from top to bot- 
tom), from 1 Mq to M,nax- Recall, that we use the model of strong 
proton superfluidity in the entire core that switches off Urea pro- 
cesses in all stars. With this assumption, neutron superfluidity (a) 
allows us to explain almost all sources, except for warmest and 
coolest ones (for their ages). However, neutron superfluidity (c) 
produces too strong CPF neutrino emission in low-mass stars and 
cannot explain the majority of warmer stars, although it agrees with 
the Cas A NS data. 

Furthermore, we can rise the cooling curves of low-mass, 
not too old NSs assuming they have more heat transparent heat- 
blanketing envelopes of light elements . For example, the dashed 
curves in Fig.|2|are calculated for M = 1 Mq NS with the carbon 



© 2010 RAS, MNRAS 000,[T]l6] 



4 P. S. Shternin et al. 




Figure 2. (Color on line) Sequences of (solid) cooling curves for NSs of masses from 1 Mq to M^^-^ (through O.IMq) with strong proton superfluidity and 
moderate neutron superfluidity (a) (left) or (c) (right) in the core (q = 0.76) compared with observations of isolated NSs. Dashed lines refer to warmest stars 
of these types - 1 M© stars with the carbon surface layer of mass 10"* M©. Dash-and-dotted lines refer to coolest Mj^a^t stars without proton superfluidity in 
the inner core. 



envelope of mass 10"** M©. We can also lower the cooling curves of 
massive (M,„ax) stars: the dash-dotted lines are computed by taking 
more realistic models for proton superfluidity, with T^pip) going 
down at high densities. This opens direct Urea process in the in- 
ner core and gives coldest possible NSs. Therefore, we can really 
explain all the data with superfluidity (a), but not with superfluid- 
ity (c). 

The second temperature drop is known in th e NS cooling the- 
ory (e.g.. iKaminker. Haensel & Yakovlevll200lh . C ooling models 
like (a ) and (c) in Fig. |2] have been analysed by lOusakov et al.l 
with the same conclusion that models like (a) can explain 
observations of all isolated NSs. Similar models of NSs with nu- 
cleon cores, where direct Urea is forbidden but C PF operates, were 
used as the basis of the minimal cooling theory jPage et al]|2004l 
I2006ll2009h although that theory employs selected T^ip) profiles, 
most favorable b y the theory of nucleon superfluidity. Now we see 
that the model of iGusakov et alj j2004) is also suitable to explain 
the Cas A NS data. 

Finally, Fig. [3] demonstrates the elfect of suppression of CPF 
neutrino emission in the axial-vector channel. To maximize the 
CPF emission we take constant over the core. It produces espe- 
cially strong splash of CPF neutrinos when the second temperature 
drop starts. Three solid lines are the cooling curves for the 1.65 M© 
star calculated at q=OJ, 0.4 and 0.19. Otherwise the conditions are 
the same as in the right panel of Fig.[T] With our constant Tc, CPF 
neutrino emission alq = 0.7 is too strong; it gives faster Cas A NS 
cooling than required by observations. The case q = 0.4 now agrees 
with the observations. Smaller q = 0.19 gives slower cooling that 
cannot explain the data. Notice that the r;.n=const model is hardly 
realistic. For more realistic T^-^ip) profiles we can reconcile theory 
with the Cas A NS data at ^ > 0.4. This gives a useful restriction 
on the uncertain theoretical parameter q. 

In the inset of Fig. [3] we show the same three cooling curves 
over larger range of ages. In addition, we plot the same dashed line 
for non-superfluid star as in the right panel of Fig. [T] and another 
dot-and-dashed line for the star with neutron superfluidity alone 
with Ten = 4.3 X 10** K and ^ = 0.76. The latter superfluidity trig- 
gers a splash of CPF neutrinos, but the main modified Urea neu- 
trino emission is too strong and the splash cannot produce a steep 
Ts decline required by the observations. Adding a carbon surface 
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Figure 3. (Color on line) Same as on the right panel of Fig.[T]but for con- 
stant Ten 

over the core at thi'ee values = 0.19 (Tc„ = 7.55 X 10* K), 
0.4 (7.2 X 10* K) and 0.7 (7 X 10* K). The inset shows the same cooling 
curves but over larger range of ages, together with the dashed curve for 
non-superfiuid star and the dash-and-dot curve for the star without proton 
superfluidity but with neutron superfluidity at T^n = 4.3 X 10* K. 



layer of mass ~ 10"'^ M©, we could raise the latter curve to the Cas 
A NS level but would be unable to reproduce the cooling slope. 
Our calculations show that the modified Urea emission should be 
suppressed at least by a factor of 30 (for the most efficient CPF 
emission with q = 0.76 and constant T^n) to get the required slope. 
Taking smaller q or narrower Ten (p) -profile would require stronger 
suppression of the modified Urea process. 



4 CONCLUSIONS 

We report a new (November, 2010) Chandr a observation of the 
young Cas A NS that confirms the observed jHeinke & Hdl201ol) 
steady decline of the surface temperature T^ (by 4% over 10 years). 
We propose a natural explanation of the observed decline. We as- 
sume that the Cas A NS underwent the traditional internal crust- 
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core relaxation some time ago and now demonstrates the second 
temperature drop due to the onset of triplet-state neutron superflu- 
idity in its core and associated neutrino emission. We can explain 
the Cas A NS observations under the following conditions: 

• The maximum critical temperature for triplet-state pairing of 
neutrons should be T^^^.^^ a (7 - 9) x 10** K. Otherwise the second 
temperature drop occurs earlier or later than required by observa- 
tions. 

• The Tc„(p) profile over the NS core should be rather wide for 
the CPF neutrino emission to gain enough strength. 

• For the same reason the suppression of the CPF process by 
collective effects cannot be too strong (q > 0.4). 

• The neutrino emission of the star before the second temper- 
ature drop should be 30-100 times lower than due to the mod- 
ified Urea process (e.g., the modified Urea can be suppressed 
by strong proton superfluidity). Otherwise the second temperature 
drop would not be pronounced. 

When these criteria are met, we can still locate T^nip) profiles 
in different parts of the NS core. If, however, we wish to explain all 
current observations of isolated NSs with one and the same T^^ip)- 
profile, we will be forced to push this profile deeper in the core 
(Fig. [2}. Alternatively, we could employ broader profiles but with 
density-dependent factor q ( which can increase within the core a s 
for singlet-state pairing, e.g.. lKolomeitsev & Voskxesenskvll2008h . 
This would shift the efficiency of the CPF process to higher p in 
superfluid matter. 

We have taken one EOS and focused on 1.65 Mq neutron star 
model but our basic conclusions will not change for a large variety 
of EOSs and masses M. For instance, taking the same EOS we have 
considered the Cas A NS models with M from 1.4 Mq to 1.9 Mq. 
By slightly changing T^„{p) profiles we are able to explain the data 
for any M from this range. 

Our calculations indicate that the second temperature drop 
lasts for a few tens of years and Cas A NS is at its active CPF 
neutrino emission stage. These models would be inconsistent with 
a sharp stop of the temperature decline in a few years, which can 
be verified with future observations. 

After the second temperature drop the Cas A NS is expected 
to become a rather cold slowly cooling NS. 
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